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1. Introduction 

An alternative title for this talk is 

What is Known About Low Energy Constants and Quark Mass Dependence in Chiral Pertur- 
bation Theory /rom the Continuum 

In order to discuss this a short introduction to Chiral perturbation Theory (ChPT) is given with an 
emphasis on some of the aspects that have been known to create confusion, the choice of quantities 
to express the expansion in. The remainder can be split in an overview of the existing calculations 
in two and three-flavour ChPT and quark mass dependences of several of the quantities of interest 
in both cases. I restrict myself here to the cases where a full two-loop calculation is available and to 
quantities of interest to the lattice community. I also remind you of the existing partially quenched 
calculations at two-loop order. In the remainder I will use order in the chiral counting, next-to- 
next-to-leading order (NNLO) and two-loop order as synonyms even though strictly speaking the 
two-loop diagrams are only part of the NNLO result. 

A review of order ChPT is ^ and several more references to lectures as well as files 
containing the long two-loop expressions can be found in the website [^. I want to point out those 
by Sharpe aimed at lattice QCD practitioners [^]. 

2. Chiral Perturbation Theory 

Chiral perturbation Theory in its modem form was introduced by Gasser, Leutwyler and Wein- 
berg [Q, ^ and it can defined as 

Exploring the consequences of the chiral symmetry of QCD and its spontaneous breaking 
using effective field theory techniques 

The assumptions that are needed to allow a derivation from QCD and a full derivation can be found 
in the paper by Leutwyler [^. I do not intend to give a full derivation here but only restrict myself 
to a few comments. 

Powercounting: In effective field theory, one assumes that there is a gap in the spectrum which 
allows to include only the degrees of freedom that are relevant below the gap and treat the effects 
of the degrees of freedom at higher scales perturbatively. Thus a clear separation of scales is the 
first requirement. One then constructs the most general local Lagrangian with the lower degrees 
of freedom in agreement with the symmetries and their realization. Unfortunately, this leads to an 
infinite number of parameters and hence there is no predictivity left. If there exists a way to organize 
the series in terms of order of importance, then we can work order by order and have predictivity 
at any fixed order in the importance. This is typically achieved by introducing powercounting. 

Unitarity: In Ref. |^ one uses strongly the fact that the only singularities at low energies in 
Green functions come from poles and cuts of the light degrees of freedom and that all the remaining 
vertices can be expanded. This is where unitarity plays a role in the derivation of ChPT from QCD. 

ChPT: ChPT is the effective field theory build with the Goldstone Bosons resulting from the 
spontaneous breakdown of chiral symmetry in QCD as primary degrees of freedom. The power- 
counting used^ is basically dimensional counting in momenta and (meson) masses. The expected 

'There are different countings possible, tlie one liere is tlie standard one. References can be found in |fn| and the 
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breakdown scale is the scale of resonances, so for energies around nip standard mesonic ChPT 
definitely does not work. The breakdown scale depnds a bit on the channel. 

Chiral Symmetry: In QCD with three light quarks of equal mass, they are fully interchangeable 
and we have a SU (3)y symmetry. But looking at the QCD Lagrangian 



q=u,d,s 



(2.1) 



we see that for = we have a separate interchange for the left and right-handed quarks: G^^ = 
SU{3)lxSU{3)r. 

Chiral Symmetry Breaking: The fact that the vacuum expectation value {qq) = {qLqR + qRqt) / 
leads to the spontaneous breaking of SU (3)l x SU (3)r to SU {3)v- The eight broken generators 
lead to eight massless degrees of freedom and makes their interactions vanish at zero momentum. 
The latter fact allows to produce a consistent powercounting in ChPT. This is illustrated in Fig. |l|. 



(;,2)2 (1/^2)2^4 





l/p' 



Jd'p 



Figure 1: An illustration of the power-counting in ChPT. On the left we see the lowest order vertex with 
two powers of momenta or masses, the meson propagator with two inverse powers and the loop integration 
leading to four powers. On the right hand-side we see two one-loop contributions and how the counting on 
the left leads to the same power p"^ for both diagrams. This counting was generalized to all orders in [Q]. 



Chiral Perturbation Theones : ChPT is a very large subject, more than 3500 papers cite at 
least one of the three basic papers. The name ChPT also is given to a variety of different theories 
exploiting the chiral symmetry of QCD. These include application to baryons, mesons and baryons 
containing heavy quarks, vector mesons and other resonances and of course the light pseudoscalar 
mesons. Within the latter we can distinguish between theories with two, three or more flavours, also 
in the partially quenched varieties, as well as including electromagnetism and the weak interactions 
nonleptonically and the possibility of treating the kaon as a heavy particle. This talk restricts itself 
to the standard two, three or more flavour sector with strong interactions and couplings to external 
currents. This is the part that has been most fully pushed to order p^. 

Lowest order Lagrangian: The Goldstone Bosons live on the manifold G/H with G = SU {np)i^ x 
SU {nF)R and H = SU {nf)v for np flavours of quarks. G/H has again the structure of SU {np) and 
can be parametrized by a matrix 



=exp(/V2<I>/Fo), with <I>(x) 



71" 



V 



K 



?r° I 178 



^/6 / 



(2.2) 



talk by S. Descotes-Genonnq]. 
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2 flavour 3 flavour 3+3 PQChPT 

F,B 2 Fo,Bq 2 Fq,Bq 2 

/ 7+3 L':,H[ 10+2 L';,H[ 11+2 

/ c't 52+4 Cf 90+4 112+3 



Table 1: The number of parameters+contact terms for the various types of ChPT. 



Here I have indicated the traceless matrix <I> in terms of the known pseudoscalars for the case of 
np = 3. The lowest order or order Lagrangian contains two terms 

■^2 = { (D^U^D^'U) + iz'^U + xU'^) } , (2.3) 

with a covariant derivative D^U = d^U — ir^U + iUl^ , which includes the left and right external 
currents: r{l)j^ = + {—)a^. The scalar and pseudoscalar external densities are included via 

X = 2Bq{s + ip) with the quark masses included in the scalar density: s = ^ H . The notation 

(A) is the trace over flavours Trj; (A). 

NLO Lagrangian: The Lagrangian at order was classified in [||, |^ and is for np = 3 

^4 = Li{D^U^D^'Uf + L2{D^U^DyU){D^'U'^D''U) +L^{D^'U'^D^UD''U^DvU) 

+L4{D^U^Df,U) ix'^U + XU^) +L5{D^U''D^Uix^U + U^)) +U{x''U + xU^? 

+Lj{x^U - xU^? + U{x'^Ux^U + xUhU^) - ih){F^yD^UD''U'' +Fl;^D^U^D''U) 

+L,o{U'^F^,UF'^n +H,{Fl^,F''^' + F;^,F'^n +H2{X^Z) ■ (2-4) 

The constants L,- in this Lagrangian are generally known as Low-energy constants (LECs). The 
constants Hi have values dependent on the definition of currents/densities and the terms are called 
contact terms. The LECs absorb the divergences of loop diagrams order by order in the powercount- 
ing. The finite part is denoted by LJ' and depends on the subtraction scale jJ. and the renormalization 
prescription. 

Number of parameters: The number of parameters at the various orders is shown in Tab. 
The order p^ is from order from ^ order p^ from [ ]T0| ] after an earlier attempt [11]. The 
partially quenched results are derived from the np flavour case [12, 13]. The difficulty in obtaining 
a minimal set can be seen from the recent discovery of a new relation in the np =2 case [Q]. 
Since the normal case is a continuous limit of the partially quenched case, the resulting LECs are 
just linear combinations of partially quenched LECs using the Cayley-Hamilton relations given in 
[jl0|]. The general divergence structure at this order is also known [15]. The parameters B^Bq and 
F ^ Fq are the two versus three-flavour lowest order constants. 

Chiral Logarithms: The main predictions of ChPT are twofold. 1) It relates processes with 
different numbers of pseudoscalars. 2) It predicts nonanalytic dependences at higher orders, often 
referred to generically as Chiral Log(arithm)s. As an example, the pion mass for = 2 is given 
at NLO by [§ 



nti = 2Bm + 



'■Tt 



1 , {2Bm) ; 



+ ••• (2.5) 



The notation = 2Bm = B{mu + md) is used a lot in the remainder. The /i dependence cancels 
between the explicit dependence in the logarithm and the implicit dependence in ly 
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LECs and choice of pL: The LECs, like Z3 in (2.5), have to be determined experimentally or 



from lattice calculations. They can be quoted in several ways. For np = 1 Ref. introduced 

/• = (327rVT4) 4'(At) - log (M^/At") , (2.6) 

which are /i independent and are proportional to the LECs Z[(/x = mj^). For np = 3 some of the 
corresponding yi are zero and no good equivalent definition of L, exists. Here we always quote the 
L\{ijl). The scale }JL is in principle arbitrary but becomes relevant when using estimates for higher 
order constants. Choosing fL = mji,mK or m-q puts some of the chiral logs to zero and thus obscures 
one of the main predictions of ChPT. At a scale /X 1 GeV experimentally Lj and this would 
clash with large Nc type estimates of the LECs. For these reasons, many ChPT practitioners use 
pL=mp= 0.77 GeV. 

What quantities to expand in: The ChPT expansion is in momenta and masses. However, one 
first has to decide whether to expand in lowest order quantities, like F, 2Brh, or physical masses 
and decay constants, like mji,mK,mr^,Fji,FK. The latter is also not unique since relations like the 
Gell-Mann-Okubo relation and kinematical relations like s + t + u = 2m| + 2mj. for Triir-scattering 
can be (and are heavily) used to rewrite expressions. This sounds trivial but can change much 
how a series convergence looks as shown below for a simple example. Similar questions are also 
discussed by [^]. I personally prefer to use physical masses and decay constants rather than the 
lowest order quantities. The physical quantities are typically better known and the chiral logs are 
created by particles propagating with their physical momentum. Also, thresholds appear in the 
right places at each order in perturbation theory. 

A very simple example: Take the relations nij^ = mo/(l +am{)/ f^) , A = /o/(l +^'«o//o) , 
as exact. We can expand to NNLO in several ways 



m^ = mo-a^ + a'^ + --- f^ = f^, I I - + b'^ + ■ ■ ■ \ (2.7) 
/o /o V /() /o / 

2 3 / 2 \ 

m„ = mo-a'^+aib-a)'^ + --- f^ = fJ I - b'^ + b{2b - a)'^ + ■ ■ ■] (2.8) 

JTt J 71 \ JTt Jn J 



m^ = mo\\-a"-^ + ab'^ + ---] (2.9) 

The coefficients in the expansion and the actual numerical values clearly depend on the way we 
write the results. The plots in Fig. ^ show the convergence for a = \, b = 0.5 and /o = 1. Only 
knowing the first three terms one would draw very different conclusions on the quality of the 
convergence from Fig. ^ for the different ways of writing the expansion. 

3. Two-flavour ChPT at NNLO 

References to order p^ and p'^ work can be found in [|l|]. The first work at NNLO used disper- 
sive methods to obtain the nonanalytic dependence on kinematical quantities, q^,s,t,u at NNLO. 



This was done for the vector (electromagnetic) and scalar formfactor of the pion in [16] (numeri 



cally) and mJU (analytically) and for ;r;r-scattering analytically in [18]. The work of mm allowed 



to put many of the full NNLO ChPT calculations in two-flavour ChPT in a simple analytical form. 
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Essentially all processes of interest are calculated to NNLO fully in ChPT starting with 77 — > 
7r°7r° [ll9|, Hi, 77 ^ n+n- ||, H, F,, and m,, [Hig HH), TTTi-scattering [||, |5]|, the 



pion scalar and vector formfactors [ ]26| ] and pion radiative decay n 
also known at order in finite volume [28]. 



Ivy ||27|]. The pion mass is 



The LECs have been fitted in several processes. I4 from fitting to the pion scalar radius [ ]27| , |29| ], 
from an estimate of the pion mass dependence on the quark masses and l\ , I2 from the 

and fg 



agreement with TTTT-scattering [|29|], Zg from the pion charge radius 
formfactor in tt ^ (vj. The final best values are [26, ^ 



from the axial 



h 

1 4 



-0.4 ±0.6, 
: 4.4 ±0.2, 



= 4.3±0.1, 
/"s = 3.0 ±0.3, 



/"s =2.9 ±2.4, 

/"e = 16.0 ±0.5 ±0.7. 



(3.1) 



It should be noticed that we do not have a good determination of It, from the continuum. 

There is also information on some combinations of LECs. These are basically via the 



curvature in the vector and scalar formfactor of the pion [ |26| ] and two combinations from TTTT- 
scattering [ ^ ] from the knowledge of bs and b(, in that reference. The order p^ LECs are 
estimated to have a small effect for mjt,fjt and TTTT-scattering. 

Let me now show a few dependences on the quark mass via = 2Bm. First for m\ expanded 



in analogy with (|2.9D. A surprise is how small the NLO and NNLO corrections are for the values of 
the input parameters given in ( ^ ) and cj'(/i = 0.77 GeV) = 0. The full result is extremely linear as 
can be seen in the left plot in Fig. ^The linearity is a consequence of the fitting parameters as can 
be seen in the right figure in Fig. ^. Similarly, Fji as a function of expanded as in (2.8 ) is shown 
in Fig. ^. The values of m\, Fji and are determined self consistently from the ChPT formulas 
quoted in [E3] via an iterative method. 



4. Three-flavour ChPT at NNLO 
4.1 Calculations 

In this section I will discuss several results at NNLO in mesonic three-flavour ChPT In general 
the formulas here are much more involved than in two-flavour ChPT and while the expressions have 
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Figure 3: The pion mass squared as a function of the quark mass viaM^ = 2Bm, left with inputs as in (3.1 ) 
and right with lo, ~ 0, both are for np —2 ChPT. 



NLO 
NNLO 



Figure 4: The pion decay constant as a function of the quark mass via m2 = 2Bm, for nf = 2 ChPT. 



been reduced to a series of well-defined two-loop integrals, the latter are evaluated numerically. 



Most recent calculations use the subtraction scheme specified in [|1^, ^ but many do not, also the 
reduction to numerical integrals is done differently by different groups which makes comparisons 



very difficult in the case of disagreement, see e.g. the discussion in |3C] about the numerical 
discrepancy with pT|]. 

The vector two-point functions were among the first calculated in the n and rj flavour quantum 



number channel [ ]33| , p2| ] and in the K flavour quantum number channel Q33I, [34j]. The isospin break- 



ing for the pco channel was done in [p5|]. The disconnected scalar two-point function relevant for 



bounds on and Lg was worked out in | |36| ] The remaining scalar two-point functions are known 
but unpublished [07|]. Masses and decay constants as well as axial- vector two-point functions were 
the first calculations which required full two-loop integrals, done in the n and Tj [|3^, ^] and the K 
channel [p 



]. The full isospin breaking contributions to masses and decay constants are in [|39|]. 
At this level many expressions were known but a full fit of LECs to experimental data could 



be carried out only after had also been evaluated to NNLO [40, 41 1. The vacuum expectation 
values in the isospin limit were done in pT]], with isospin breaking in [39] and at finite volume 
in [HI. 
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The vector (electromagnetic) formf actors for pions and kaons were calculated in [43, M, 45] 
and in [ p5| ] a NNLO fit for Lg was performed. L\q can be had from the axial formfactor in 71,^ — > 
Ivy. The NNLO calculation is done, but no data fitting was performed]^]. A rather important 
calculation is the formfactor. This calculation was done by [|^, 31 1 and a rather interesting 
relation between the value at zero, the slope and the curvature for the scalar formfactor obtained 



[30]. I will present some results for below. 



Calculations for scalar formfactors including sigma terms and scalar radii [ ]47[ ] and mz [ |48[ ] 
and TiA'-scattering [ p9| ] have been performed as well and used to place limits on L\ and Lg. Finally, 
the relations between the /[ and LJ" has been extended to the accuracy needed to compare order 
results in two and three-flavour calculations [^] and there has been some progress towards fully 
analytical results for [55] and Ti^-scattering lengths [56]. 

I am aware of number of calculations in progress, including 77 ^ 3?! and isospin breaking in 



Kb fl51|], preliminary results were reported in [|52|], the sunsetintegrals needed for the masses at 
finite volume [|53|] and the relations between the two and three flavour order constants [54]. 



4.2 C[: estimates of order p^ LECs 

Most numerical analysis at order p^ use a (single) resonance approximation to the order p^ 
LECs. This is schematically shown in Fig. ||. The main underlying motivation is the large Nc limit 



p,5 



71 



71 




Figure 5: A schematic indication of the estimate of the order p LECs by resonance exchange. 



and phenomenological success at order p'^ [57|. There is a large volume of work on this, some 
references are [58, 59, 6C, ^]. The numerical work I will report has used the simple resonance 



Lagrangian [ |25| , [39[ , |41 



1 



2V2 



1 



1 



= -{y''SV^S-Mis')+Cd{Su^'u^)+c,„{Sx+), ^r^f = -d^P^d^Pi- -M'^,Pt + idMx 



/y = 0.20, /;^ = -0.025, gv=0.09, c„, =42MeV, = 32 MeV, 4 = 20 MeV, 



my 



■■ 0.77 GeV, Ma = rua, = 1 .23 GeV, ms = 0.98 GeV, mp^ = 0.958 GeV . 



(4.1) 



The values of /y, gv, fx and come from experiment [ |25| , p7p and and from resonance 
saturation at order p"^ [^]. 

The estimates of the is the weakest point in the numerical fitting at present, however, many 
results are not very sensitive to this. The main problem is how the C[ which contributes to the 
masses, estimated to be zero except for r]' effects by (pTT|), affect the determination of the others. 
The estimate is also /i -independent while the C[ depend on }x. 
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The fits done here in [39, 40, 41, ^] try to check this by varying the total resonance con- 
tribution by a factor of two, varying the scale /i from 550 to 1000 MeV and compare estimated 
C[ to experimentally determined ones. The latter works well, but again the experimentally well 
determined ones are those with dependence on kinematic variables only, not ones relevant for 
quark-mass dependence. 



4.3 The fitting and results 



The inputs used for the fitting, see the more extensive discussion in [39, 41], are 



Ki^\ F{0), G(0), A from E865 at BNL[|2[ |3|]. 

m^o, m^, m^o> electromagnetic corrections include the estimated violation of Dashen's 
theorem ([64] and references therein). 



nis/m = 24. Variations with m,f/rh were studied in [||,||]. 

L4,Lg the main fit, 10, has them equal to zero, but see below and the arguments in [36] 



The results of this fit are summarized in Tab. ^ The errors are very correlated, this is shown in 
Fig. 6 in [^] for the fit to then available Kf4 data. As said before varying the estimates of the C[ by 
a factor of two, varying the /i where the estimate of the is done from 550 to 1000 MeV all stay 
within the given errors. Varying the values of L4,Lg as input can be done with a reasonable fitting 
chi-squared when varying lO^L^ from —0.4 to 0.6 and from —0.3 to 0.6. These alternative fits 



were performed in [47] and variation of many quantities with L4,Lg (including the changes via the 



changed values of the other LJ") are shown in [47, ttq, 49]. Fit B was one of the fits that gave a good 



fit to the pion scalar radius and fairly small corrections to the sigma terms [ |47| ] while fit D [ ]650 is 
the one that gave agreement with nn and Tr/T-scattering threshold quantities. 

One point should be observed here, if one fits lattice data with NLO formulas to obtain the L-", 
one should also use the NLO or order fit values from Tab. ^ to compare with. 

Note the niu/m^i = is never even close to the best fit and this remains true for the entire 
variation withL4,L^. The value of Fq, the pion decay constant in the three-flavour chiral Umit, can 
vary significantly, even though I believe that fit B is an extreme case. 

In Fig. ^ we show how the threshold parameters ajj and depend on the variation with Z.4,Lg. 
always agrees well with the result of [29] while Oq only agrees well within a limited region [p^]. 



For comparison, the order values are aQ = 0. 159 and = —0.0454. The planes in Fig. ^ indicate 



the results 



0.220 ±0.005, al = -0.0444 ±0.0010 [g9p. The same study was performed for 
TiK scattering lengths in [ p9| ] with the results of the Roy-Steiner analysis [^6|]. The resulting limits 
on the input values of L'^,L^(, are shown in Fig. 0. The resulting region called fit D in Tab. § is 
10^ L4 « 0.2, lO^L^ 0.0. This general fitting obviously needs more work and systematic studies 
and constraints from lattice QCD on L4,Lg will be very useful. 



4.4 Mass dependence and other results for selected quantities 

I now show the dependence of a few quantities on the input masses. These are updates of 
the plots shown in [ pl| ] and also some new ones for /+(0) in Ki^. The masses squared and decay 
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fit 10 


same p^ 


fitB 


fitD 




0.43 ±0.12 


0.38 


0.44 


0.44 




0.73±0.12 


1.59 


0.60 


0.69 




-2.53 ±0.37 


-2.91 


-2.31 


-2.33 




= 


= 


= 0.5 


= 0.2 




0.97 ±0.11 


1.46 


0.82 


0.88 




= 


= 


= 0.1 


= 




-0.31 ±0.14 


-0.49 


-0.26 


-0.28 




0.60±0.18 


1.00 


0.50 


0.54 




5.93 ±0.43 


7.0 


_ 


— 


2Bom/m^ 


0.736 


0.991 


1.129 


0.958 


2 4 6 

P ,P 


0.006,0.258 


0.009,= 


-0.138,0.009 


-0.091,0.133 


2 4 6 


0.007,0.306 


0.075,= 


-0.149,0.094 


-0.096,0.201 


2 4 6 


-0.052,0.318 


0.013,= 


-0.197,0.073 


-0.151,0.197 


rriu/md 


0.45±0.05 


0.52 


0.52 


0.50 


Fq [MeV] 


87.7 


81.1 


70.4 


80.4 




0.169,0.051 


0.22,= 


0.153,0.067 


0.159,0.061 



Table 2: The fits of the LJ and some results, see text for a detailed description. They are all quoted at 
/i = 0.77 GeV. Table with values from [|9[ |5|, |7[|^ HJ. At present the best fits to use for compai'ison 
with the lattice are fit 10 at NNLO or order p'^ depending on the order of the lattice fit. 



constants are written in the form analogous to (2^) as published in [33] (note the erratum of [39] 
and the formulas given in [^). A selfconsistent set of wi^, m^, m^,F^, Bonis and Bom with the fitted 
values of L[ and Fq is determined for each input value of two masses. This is done by iterating the 
formulas till convergence is reached. NNLO reproduces the physical values at the physical point. 

We show for fit 10 and fit D keeping nis/m = 24 and varying in Fig. |8[ The same 
dependence but for mj^ is shown in Fig. ^. The large corrections for fit 10 come from the kaon 
mass. This is shown in Fig. [T^a) where we plot ni^ with m fixed to its physical value and vary nig. 

The decay constants and ratios are shown as a function of nis at fixed ratio rris/m for F„ in 
Fig. [igCb), Fk in Fig. |ll](a) and Fk/F^ in Fig. |lT](b) 



Let me finally discuss some results of the calculation of Ka [l30|]. One major point is that the 
scalar formfactor can be written as 

/o(0 = 1 - 4 {C[2 + K - ml) ' + ^ {2C[2 + C,,) {ml + ml) - ^t^C'u + A(0 .(4.2) 
^jt ^7c 

A{t) contains no C\ and only depends on the at order p^. All needed parameters C\ for f+{0) = 
/o(0) can thus be determined experimentally or from the lattice via the slope and curvature of /o(f). 

The dependence of /+ (0) — 1 on the masses is shown in Fig. [l^. The left plot shows the three 
corrections, order p'^, order p^ pure two-loop part and order p^ the LJ"-dependent part as well as the 
sum of those three. It does not include the contribution from the C\. The right plot shows the same 
contributions divided by {m\ — Both ai^e shown as a function of m| with m\ fixed. Note that 
in the right plot both order p^ corrections are rather fiat except near the physical pion mass. 
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Figure 6: The nn scattering lengths in three-flavour ChPT as a function of the input values of ^4,^5 used in 
the fits, fljj left, flg right. See [Q for details, from [Q. 



jTjt constraints tiK constraints 




'0.4 -0.2 0.2 0.4 0.6 -0.4 -0.2 0.2 0.4 0.6 

10^ l; tO^i 

Figure 7: The bounds on L^^U^ from nn and ;r/r-scattering threshold parameters. Left nn where the bound 
from Oq shown in Fig. ^ is the most stringent. Right TzK. White regions are allowed. The region of fit D, 
compatible with both, is indicated by the circle. From [|4^]. 

5. Even more flavours at NNLO (or PQChPT) 



NLO Partially Quenched ChPT has been studied by many people and found to be very useful, 
see ^ and references therein. The masses and decay constants are known to NNLO for almost 
all possible mass combinations. Formulas were kept in terms of the quark-mass expansion, analo- 
gously to (17), to avoid the proliferation in physical masses appearing in this case. The three sea 
flavour masses and decay constants can be found in [12, 13, 57] and the two sea flavour results are 
in [58]. Numerical programs are available from the authors. The formulas are in the papers but can 
be downloaded from [^. The papers also contain discussions about how to fit the NNLO LECs. 
PQChPT NNLO results for neutral masses are in [E^] and electromagnetism was included in [7C]. 
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Figure 8: m\ as a function of for fit 10 (left) and fit D (right) of Tab. g with nis/m fixed. Note the 
difference in convergence properties between the two fits. 





Figure 9: m\ as a function of lUs for fit 10 (left) and fit D (right) of Tab. || with nis/m fixed. Note the 
difference in convergence properties between the two fits. 




(a) 



NLO ■ 
NNLO - 



(b) 



Figure 10: For fit 10: (a) as a function of with m fixed, (b) F-,t as a function of nis with nis/m fixed. 
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6. Wishing lists and Conclusions 

The conclusions are several lists. First a general wishing list and a comment on NNLO fitting: 

• Quark mass dependences everywhere 

• A presentation of lattice results at a given quark mass extrapolated to the continuum and 
infinite volume. Or, in general, results presented in a way that allows us ChPT practitioners 
to use your data also later on when other inputs might be changed. 

• More use of the existing NNLO calculations, if you have any ideas how we can improve the 
usability of the existing calculations, please tell us. But remember: 

- The number of new parameters at NNLO is very small compared to present lattice 
fitting results. The new order parameters are exactly the same in number as those 
you add when you add an analytic NNLO fitting expression. You therefore do not loose 
any predictivity when including the full NNLO result and including the known NNLO 
parts is definitely recommended. 

- There is a small caveat to this, the order LECs that show up in scattering now appear 
at NNLO also for the masses and decay constants but these LECs are well known for 
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Hf = 2, l\,l2 from Eq. (3.1), and np = 3, L\,L2,L'^ from Tab. |^ For the partially 



quenched case, there is the unknown extra parameter Lq, but a large Nc estimate gives 



L[=L'2 = and thus = 1L\, L\ = Ui+'^L\, (obtained by inverting (21) of [|1J). 



A more direct wishing list (with input from Berne) for two-flavours is 

• Z3 and errors. 

• l^: can the lattice check the relation between Fji as a function of the masses and the scalar 
radius? 

• U and [29 1 Can the lattice check the strong correlation between the scalar radius and a^l 



• Isospin breaking in kk scattering at 5 = m^, important for CP-violation phenomenology. 

A similar list for three flavours (the different points are rather related) 

• Large N^- suppressed couplings like L4 and 

• wis dependence of Fji and m\ 

• sigma terms and scalar radii 

• iri J^B from extrapolations of both f+{q^) and fo{q^) 

• sigma terms and scalar radii 

and my final comments 

• Lots of analytical work is done in ChPT at NNLO, please use (and cite ©) it. 

• Use the correct ChPT 

- 2-flavour for varying th and possible for Nf = 2 and Nf = 2+ I at fixed (but have 
different LLCs) 

- otherwise 3-flavour 

- the various partially quenched versions 

• Remember at which order in ChPT you compare things: NLO for both lattice and continuum, 
NNLO for both lattice and continuum. Fits put neglected higher order effects into the LECs. 

And finally, ChPT and LECs played to my great pleasure a large role in many of the presenta- 
tions at this conference, as summarized by S. Neccoff/ll]. I am looking forward to even more future 
results from lattice QCD. 
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